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Abstract. We report on the generation of non separable beams produced via the 
interaction of a linearly polarized beam with a cloud of cold cesium atoms placed in 
an optical cavity. We convert the squeezing of the two linear polarization modes into 
^vq ■ quadrature entanglement and show how to find out the best entanglement generated 

in a two-mode system using the inseparability criterion for continuous variable [Duan 
et ai, Phys. Rev. Lett. 84, 2722 (2000)]. We verify this method experimentally with 
a direct measurement of the inseparability using two homodyne detections. We then 
map this entanglement into a polarization basis and achieve polarization entanglement. 
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PACS numbers: 42.50.Lc, 42.65.Pc, 42.50.Dv 
1. Introduction 

With the recent progress in the quantum information field, there has been a lot of 
interest in entanglement in the continuous variable (CV) regime. Criteria to demonstrate 
and quantify CV entanglement have been developed [JJ Ej and experimentally tested 
[3 IU El HI IE] In particular, the possibility to map a quantum polarization state 
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of light onto an atomic ensemble j2] has stirred a great deal of attention to the 
quantum features of polarized bright beams. The notion of polarization entanglement, 
i.e. entanglement between Stokes parameters of two spatially separated beams, has 
been investigated by Korolkova et al. [10J and first demonstrated by Bowen et al. 
[7j by mixing two independent squeezed beams produced by OPAs. Polarization 
entanglement was also achieved via the Kerr non linearity of optical fibers [IT] and 
cold atoms jHJ- These experiments are important steps in connection with quantum 
teleportation |T2|, quantum dense coding [13], entanglement swapping [13] and, more 
generally, characterizing entanglement in the CV regime |15j . 

In this paper, we report on the generation of non separable beams via the interaction of 
a linearly polarized light beam with a cloud of cold cesium atoms placed in an optical 
cavity [H|. In previous works [TB] ll7j. we have shown that, after the non linear interaction 
with the atoms, two modes of the light exiting the cavity were squeezed: the mean field 
mode, but also the orthogonally polarized vacuum. We develop here a general method 
to find out the best entanglement - as measured with the inseparability criterion [2] - 
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produced in a two mode system and characterize the correlation properties of the system 
in the Poincare sphere. The main result is that the maximal entanglement corresponds to 
the sum of the minimum noises of two " uncorrelated" modes. The maximally entangled 
modes are then circularly polarized with respect to these modes. We stress the similarity 
with the usual entanglement experiments, which mix independent squeezed beams on 
a beamsplitter [HE]- Moreover, we show that this mixing of two independent beams 
is equivalent to rotating the polarization basis of a single beam exhibiting correlations 
between polarization modes. This intuitive approach could be of interest for the study 
of systems in which quantum correlations exist between polarization modes, and allows 
one to think in terms of independent beams mixing. 

We then apply these results to our experiments [Sj, show experimental evidence of 
both quadrature entanglement and polarization entanglement. In Sec. El quadrature 
entanglement is demonstrated by figuring out the maximally entangled modes and 
checking the inseparability criterion between these modes in a direct detection scheme. 
In Sec. E] we map the entanglement into a polarization basis via the mixing of our 
quadrature entangled modes with an intense coherent beam, the phase of which is 
locked to that of the first beam. We therefore measure directly the Stokes parameters 
fluctuations of the two spatially separated beams, thus demonstrating polarization 
entanglement. 

2. Looking for maximal entanglement 

2.1. General method 

In this section we develop a general method to find out the maximal entanglement in 
a two-mode system. We start with a "black box" - in our case the atomic medium in 
the cavity - out of which comes a light beam with unknown quantum properties. Let 
us stress that the goal of this Section is to develop a method to characterize quantum 
properties, such as entanglement and squeezing, which have been previously created 
between some polarization modes by some interaction. Let us denote by A a and At, 
two orthogonally polarized modes of this beam. They satisfy the standard bosonic 
commutation relations [A a , AL] = 5 a p- The usual quadrature operators, with angle 9 in 
the Fresnel representation, 



are the continuous variable analogous of the EPR-type operators as introduced by 
Einstein, Podolsky and Rosen ^H]- The criterion derived by Duan et al. and Simon [2] 
sets a limit for inseparability on the sum of these EPR-type operators variances 



For Gaussian states, X a ^(9) < 2 is a sufficient condition for entanglement and 
has already been used several times to quantify continuous variable entanglement 



X a (9) = A a e- ie + Ale 



Y a {6) = X a (9 + tt/2) (a = a,6) 




(1) 
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[HH H21 HH El- In this Section we look for the best entanglement produced in the 
system: using unitary transformations we therefore seek to minimize I ajb with respect 
to a, b and 9. Expanding (JTJ), one gets 

l a , b {9) = (5A{5A a +5A a 5A{+5Al5A b +5A b 5Al)+A\(5A a 5A b )\ cos[2(0-0 a , 6 )](2) 

where 9 a ^ b is the phase of (SA a SA b ). The minimum value is reached for 9 = 9 a)b ±7i/2: 

X afi = mml a , b (9) = (6A{6A a +6A a 5A{+6Al6A b +5A b 5Al) - 4\(5A a 5A b )\(3) 

6 

X ab does not depend on local unitary operations performed separately on a, b. It 
thus provides a good measurement of the entanglement between modes a and b and will 
be used throughout this paper. Consequently, one has to look for the polarization basis 
(a*, b*) of the "maximally entangled modes" which minimizes T a b . It is easy to see that 
the first term in (jHJ) is independent of the polarization basis, since it is the trace of the 
correlation matrix of modes a and b. The entanglement between a and b is therefore 
completely determined by the correlation term | (SA a SA b ) |. 

In order to find the strongest correlations we turn to a particular basis for the 
fluctuations (we are only interested here in what happens to the "noise ellipsoid", 
regardless of the mean field): as shown in the Appendix A, there always exists two 
orthogonally polarized modes A u and A v such that (SA U 5A V ) = 0. The u, v modes are 
"uncorrected" in the sense of the inseparability criterion and satisfy X UjV = m&x ajb I atb > 
2. Note that these modes are not uncorrelated stricto sensu, since (5A u 5Al) can be non 
zero. Moreover, our choice is not unique, since any A' u = e idu A u and A' v = e V A V also 
satisfy the same property. To unambiguously determine the "uncorrelated" basis we 
choose modes u and v such that (SAD and (5A 2 ) are positive numbers. Physically, it 
means that we choose u and v such that their noise is minimum for the same quadrature 
Y. 

Two orthogonally polarized modes a and b decompose on such a basis 



A a = (3A U - ae^A v (4) 
A b = aA u + (3e^A v (5) 

with a, (3 positive real numbers such that a 2 + (3 2 = 1. The correlation term reads 

\(5A a SA b )\ 2 = a 2 p 2 {(6A 2 u ) 2 + (5A 2 v ) 2 -2(6A 2 u )(5A 2 v )cos2<i ) } (6) 

and is maximal for <\> = n/2 [tt] and a = (3 = 1/V2. The maximally entangled 
modes are then the circularly polarized modes with respect to modes u,v: 

A a * = -j= ( A u - iA v ) (7) 
A b » = -j= ( A u + iA v ) (8) 
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which satisfy 

\{SA a *5A b *)\ = m a x\(5A a 5A b )\ = ~[ (5A 2 U ) + (5A 2 V ) } 

a,b Z 

Plugging this result in Q and using (J7HHJ), we obtain the maximal entanglement 
value as the sum of the minimal noise of the "uncorrelated" modes 

Z a * )6 * = minJ aj6 = (5X 2 a ) rnin + (5X 2 ) min (9) 

a,b 

This is the main result of this Section: the best entanglement in the system is found 
between the circularly polarized modes in the (u, v) basis, and it is equal to the sum of 
the u, v mode minimal noises. This result stresses the link between entanglement and 
squeezing; associating to each polarization basis (a, b) the sum of the minimal noises 

S a ,& = (5X 2 ) min + (5X 2 ) min (10) 
= {8A\8A a + SAJAl + 8A\8A b + 5A h 5A\) - 2[ \(8A 2 a )\ + \{5A 2 )\ ] 

Eq. © clearly shows that looking for maximal entanglement is equivalent to looking 
for maximal squeezing produced by the system 

l a *,b* = min \(5X 2 ) min + (5X 2 ) min \ = T, m i n (11) 

a,b 

2.2. Correlations in the Poincare sphere 

A standard representation for the polarization state of light is provided by the Poincare 
sphere [10], which relies on the Stokes parameters [20]. Given the special role played by 
the "uncorrelated" basis, we define the Stokes parameters from the u, v modes 

a' — a* a _i_ 4*4 o'_4*4 _ A* A 

S 2 = A^Ay + A*A U 5*3 = i(A*A u — A* U A V ) 

and we study the evolution of entanglement and squeezing when the polarization 
basis is rotated, that is, when the polarization state vector moves along the Poincare 
sphere. In the general case the correlation properties of the system can be summarized 
as follows (see Appendix B for the demonstration of these results): 

(i) Along the "uncorrelated" modes axis S[, I a ,b is maximal by construction 
((5A U 5A V ) = 0) and these modes are never entangled 

I UtV = max I a:b = (5A j JA u + 5A U 5A{ + 5A\5A V + 5AJA\) 

a,b 

= (5X 2 U ) + (5Y 2 ) + (5X 2 ) + (5Y 2 ) > 2 




B u =max{fi aJ ,} 

Figure 1. Quantum properties of the beam in the Poincare sphere. The arrows with 
sign "+" correspond to increasing correlations. 

These modes are characterized by the fact that X UyV (8) is independent of 6. The 
least noisy quadratures are the same and the noise reduction is maximal 

(ii) In the equatorial plane (S[, S' 2 ), corresponding to the linearly polarized modes with 
respect to u,v, the noise reduction is also maximal and equal to Tj min . However, 
the entanglement is not constant: the best entanglement is obtained along the S' 2 
axis (modes at 45° to the u, v modes) and its value is equal to the weakest noise 
reduction H max 

1 c ,d — ^max = min [(SX 2 ) min + (SX 2 ) max , (SX 2 ) max + {SX 2 ) min \ 

(iii) Along the S' 3 axis are the maximally entangled modes a*,b*, for which the 
entanglement is maximal (and equal to the best noise reduction value) 

I'd* ,b* mifr-^a,ft ^min (^^ u )rran {^-^- v )min 
a,b 

and for which the excess noise is the largest 
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A graphical representation of these results is given in Fig. ^ at the poles there is 
maximal entanglement (T a * b * = S min ) and the worst noise reduction (£ a * j6 » = S ma x)- 
In the equatorial plane the noise reduction is optimal (£sj,s^ = S min ). Note that the 
c, d modes at "±45°" may be entangled while having optimal noise reduction. 

2.3. Interpretation 

An interpretation of the previous results can be given in connection with typical 
experiments in which entanglement is obtained by having independent beams interfere 
with each other pHE]. Indeed, going from the (u, v) to the (a, b) basis with parameters 
a, P and [Eqs. f!4|5|l ] is equivalent to combining the u, v modes on a beamsplitter with 
transmission T = /3 2 , the v mode being dephased by <p [see Fig. Efa)]. Although, again, 
the u, v modes are not completely uncorrelated, we can interpret the general results 
enunciated above as the result of an interference between two independent beams. This 
configuration is the one typically used to generate entangled beams: two squeezed beams 
on the same quadrature are produced separately - for instance with two OPAs [THE], 
or by using the Kerr-type non-linearity of optical fibers P~Ij - and then combined on a 
beamsplitter. 

The u,v modes are squeezed for the same quadrature Y, so that, in case of a zero- 
dephasing before the beamsplitter, the squeezed and noisy quadratures do not mix 
and the noise reduction sum is conserved [Figl^b)]. The outgoing beams are linearly 
polarized with respect to the incoming ones (<f) = 0) and we retrieve the property (ii) 

If the dephasing before mixing is now equal to 7r/2, the outgoing beams have excess 
noise on both quadratures [see Fig. I2fc)]. However, these noises are correlated and, 
obviously, all the more so for a 50/50 beamsplitter, since it maximally mixes the u, v 
modes. This transformation is equivalent to going from the linearly polarized basis (u, v) 
to the circularly polarized basis (a*, b*), which is naturally the correlated basis 

\{5{X a , + X„?{6)) = (5X 2 U (6)) < 1 f<* 9 = 1 

\(5{Y a , - Y h *)\6)) = (5X 2 V (6)) < 1 far = 1 
Thus T a * jfc , = (5Xl) min + (5X 2 u ) min < 2 

When the u,v modes are symmetrical (same noise properties), it can be shown 
that they are completely independent (see Sec. El for an example). The a*,b* modes 
fluctuations are then the same for all quadratures. 

Before applying these results to our experiment, we would like to emphasize that 
this analogy between the quantum properties of any system and those produced via 
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Figure 2. Interpretation of the correlations. (a) Equivalence between the 
transformations corresponding to the transmission by a beamsplitter (left) and a 
polarization basis rotation (right), (b) Squeezing conservation for the linearly polarized 
modes when there is no dcphasing between the squeezed quadratures, (c) Generation of 
entangled modes via the interference of two modes squeezed for orthogonal quadratures. 



the mixing of two independent beams provides us with a simple interpretation of the 
main results. This analysis is of particular interest to the study of systems for which 
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the correlations are simultaneously produced inside a "black box". One has to find 
the " uncorrelated" modes u, v, an experimental signature being that the entanglement 
value I u>v (9) does not depend on 9. Once this basis is obtained, one can apply the 
previous formalism or, equivalently, think in terms of independent beams combination. 
This method allows, for instance, to determine the maximal correlations produced by 
an OPO inside which are inserted birefringent elements [Hi] . 



3. Quadrature entanglement 

In our experiment J7j, the black box consists in an optical cavity containing a cloud 
of cold cesium atoms into which is sent an off-resonant light beam. In this Section we 
first present the principle of a detection scheme allowing for a direct measurement of 
the quadrature entanglement given by (JTJ). We then study the entanglement generated 
in the case of a linear incident polarization, which is qualitatively different from the 
circular polarization case. 



3.1. Entanglement measurement principle 
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Figure 3. Schematic of the a, b modes entanglement measurement. 



In order to measure the entanglement given by ([TJ between two orthogonally 
polarized modes a and b, we can reexpress X a ^(9) as the sum of the noises of circularly 
polarized modes 1, 2 with respect to a, b 



l a , b (9) = (6X1(9)) + (5X 2 2 (9)) 
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V2 



{A a + A b 



A h 



As represented in Fig. El the modes 1, 2 are straightforwardly obtained from the 
given a, b modes with a half-wave and a quarter-wave plate. They are then mixed 
with a strong coherent local oscillator (LO) on a polarizing beamsplitter and sent to 
two balanced homodyne detections. We thus simultaneously measure the spectral noise 
densities (5Xf(9)) and (<5A|(#)) at a given analysis frequency, the sum of which directly 
gives l ay b(0). This value oscillates when the LO phase is varied in time. Note that, 
unlike usual detection schemes jTj involving two successive measurements, this method 
is based on one simultaneous measurement. 

3.2. Case of a linear incident polarization 
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Figure 4. Experimental set-up. 



In the system considered in [8], an x-polarized beam interacts with a cloud of cold 
cesium atoms in an optical cavity. The experimental set-up is shown in Fig. The 
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cavity coupling mirror has a transmission coefficient of 10%, the rear mirror is highly 
reflecting. We probe the atoms with a linearly polarized laser beam detuned by about 
50 MHz in the red of the 6S1/2, F=4 to 6P3/2, F=5 transition. The optical power of 
the probe beam ranges from 5 to 15 /zW. After exiting the cavity, both the mean field 
mode A x and the orthogonally polarized vacuum mode A y are squeezed for frequencies 
ranging between 3 and 12 MHz. An interpretation of these results ^H] can be provided 
by modelling the complicated 6S1/2, F=4 to 6P 3 / 2 , F=5 transition by an X-like four-level 
atomic structure (Fig. |HJ). 




Figure 5. Atomic level structure considered: X-like configuration. 



When the two transitions are symmetrically saturated, the atoms behave as a Kerr- 
like medium for the circular components 



At = -^(A =F iA y ) (12) 

both of which are squeezed ^E]- Because of the symmetry of the system, they are 
obviously squeezed for the same quadrature. From the viewpoint of linear polarizations 
the x, y modes are also squeezed due to cross-Kerr effect, but for orthogonal quadratures 

Indeed, from (JUJ), one derives the following expressions 



(SAJAy) = -^((SA 2 + ) - (5A 2 _) + (SA+SAi) - (SA^SA-)) = (13) 
v2 

(6AJAD = ^{(5A + 8A\) - (5A_5Al) + (SA+SAl) - (5Al5A_)) =(101) 
V2 

(SAD = ^= ((SAD + (SAl ) - 2(5A + 5A_)) (15) 
(SAD = - ^((SAl) + (SA 2 _) + 2(5A + 8A.)) (16) 

Eqs. f|13H14|) show that the x, y modes are completely independent. We then 
measure X Xty (0) following the previous procedure: A +45 and iv4_ 45 are sent to the 
homodyne detections (Fig. HJ), yielding the quantity 
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Figure 6. X x ^ y {9) when is varied in time, at 5MHz. 

We verify in Fig. El that this quantity is effectively independent of 8, ensuring that 
the x,y modes are " uncorrelated" in the sense of the inseparability criterion. 
Bearing in mind that the u, v mode noises must be minimal for the same quadrature, 
we look at (|15H16|) to find the relative orientation of the minimal quadratures of the x, y 
modes. These minimal quadratures are a priori different and depend on the correlation 
term (8A+8AJ) between the circularly polarized modes. Physically, these modes are 
correlated by optical pumping processes between Zeeman sublevels. However, we place 
ourselves in the "high frequency" limit for which the analysis frequency (a few MHz) is 
much higher than the optical pumping rate (a few hundreds of kHz). The modes A± 
are then uncorrelated ((8A + 5A-) ~ 0), and squeezed for the same quadratures. Using 
(IToTfToT) and the fact that 



(8X^(9)) + (5Xf { _ m (9)) 

^(AC + Ay) 



(17) 



V2 



(A, 



A,, 



{8X1) = (SAJAl + 6AlSA a ) + 2(8 A 2 a ) cos 26, (a = u, v) 

we deduce that the x, y modes are squeezed for orthogonal quadratures. These 
properties were verified both theoretically ^B] and experimentally |T7j [see also Fig. 
|7[a)]. In the high frequency limit, one has thus to dephase one mode by n/2; we choose 
for "uncorrelated" basis A u = A x and A v = iA y . 

The maximally entangled modes are then the modes at 45° to the x, y basis 

A a * = ^=(4. - iA v ) = — {A x + A y ) = A +45 (18) 



A b . =-j=(A u + iA v ) = -j={A x - A y ) = A_ 45 (19) 
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and their entanglement is measured by summing the minimal noises of A x and iA y 




The results are reproduced in Fig. [7| for an analysis frequency of 5MHz: the 
two modes have indeed the same spectrum, and, when summed, the minimal value of 
1+45,-45(6*) is below 2, demonstrating entanglement. Equivalently, one could have set 

A u — stud -Ay — A ? since these modes are uncorrelated and symmetrical. To 

obtain the entangled modes, one has to dephase them by 7r/2 and combine them on a 
beamsplitter, yielding again the A± 45 modes 



1 e*5 



A_ 45 = —{A x - Ay) = j={A + - iA_) 



1 

Let us now represent these results in the Poincare sphere. In contrast with the 
previous Section, we define the Stokes parameters in a more usual fashion from the 
linearly polarized modes (2, y) basis 
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Figure 8. Entanglement for the canonical polarization basis: (x,y) (a), (<7+,<7_) (b), 
(+45,-45) (c), for an analysis frequency of 5MHz. 



A* A — A* A 



S = A* X A X + a;a v s 1 
1S2 = A x Ay + A A x S3 = i\AyA x A x A y ) 



With the relations A u = A x , A v = iA y , it is easy to see that the new Poincare 
sphere is obtained from the previous one via a rotation by tt/2 around the Si axis. If 
we measure the entanglement in the (cr+, a—) basis 



and represent the results in Fig. |H1 one sees that the basis in the (Si, S 3 ) plane are 
uncorrelated, as expected, whereas maximal entanglement is found for the ±45° modes. 
We also measured the noise spectra of each modes in the three basis and checked that 
the squeezing is maximal and identical for all modes in the (Si, S3) plane 



■mm 



(5X 2 a+ ) min = (5X 2 a _) min ~ 0.95 
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and that the entangled modes A± i5 have almost isotropic fluctuations in the Fresnel 
diagram, as well as identical spectra because of the independence of the x, y modes. 
One can summarize these results in the Poincare sphere represented in Fig. |H1 




Figure 9. Quantum properties of the beam in the Poincare sphere at high frequency, 
in the linear polarization case. 




Figure 10. Quantum properties of the beam in the Poincare sphere at lower frequency ; 
in the linear polarization case. 

3.2.1. Frequency dependence The situation is a little bit more complicated at lower 
frequencies. The correlations between the circularly polarized components must be 
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taken into account for frequencies lower than the optical pumping rate. The <r± modes 
are still symmetrical and the x, y modes still uncorrelated. However, the minimal noise 
quadratures are rotated [Eqs. (jl5H16j) ] and the x, y modes are no longer squeezed 
for orthogonal quadratures. To retrieve the "uncorrelated" basis, one mode must be 
dephased, say y: 



2\(5A+5A-)\mn(<l>2-<f> l ) 



where tan 0c 



(SA\) + (SAL) + 2\(6A+6AJ)\cos{<f> 1 - 2 ) 
In the high frequency limit, |((L4 + (L4_)| — > and, therefore <pc ~^ 0: we retrieve 
the previous (u, v) basis. This dephasing for the y mode is equivalent to a rotation of the 
Poincare sphere by an angle <pc around the Si axis [Fig. ^j. The maximally entangled 
modes are then 



A a * = -^={A U - iA v ) = ±={A X + e-^Ay) (21) 



A b * = —(A u + iA v ) = —{A x - e^Ay) (22) 



^ {Au + lAv) ~72 

and their entanglement is still given by (j2"Uj) 



Z a *,b* — (SX x ) min + (SX y ) min 

In Fig. [TT]are plotted the x, y modes squeezing versus frequency, as well as 1+45^45 
and the optimal entanglement T a *^*- At low frequency the squeezing improves for the 
vacuum mode A y , but degrades for the mean field mode x, so that the entanglement 
actually decreases at low frequency. In Fig. [TT]we also report the value of X +45j _45. One 
sees that it equals the optimal entanglement in the high frequency limit, but, for lower 
frequencies, the two values differ, confirming that the maximally entangled modes are 
no longer v4 ±45 , but given by f)21M22|) . 

3.3. Case of a circular incident polarization 

In this Section, we illustrate the differences between a two squeezed mode system and 
a single squeezed mode system. We show that the entanglement produced in the 
latter is qualitatively different from the former, even if the amount of correlations is 
the same. Whereas the two squeezed modes situation corresponds to the case of a 
linear polarization, the single squeezed mode situation appears when the polarization 
is circular. Indeed, our system may exhibit polarization switching: the intracavity 
polarization may become circular under some conditions [Tfi] . In this case, the atoms 
only interact with one mode - say cr + , which may also be squeezed because of Kerr effect. 
Yet, the situation is very different from the previous one. One can set A u = A + and 
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Figure 11. (a) x, y modes minimum noises versus frequency, (b) Entanglement for 
the ±45°-polarized modes (T + ^ 5 ^_^) compared to the maximal entanglement (1 a *,b*)- 



A v = A_, the former being squeezed, the latter being a coherent vacuum. The vacuum 
fluctuations are isotropic and, therefore, the properties of the (w, v ) basis must remain 
unchanged when the v mode is dephased. In other words, the Poincare sphere must be 
invariant under rotations around Si. Moreover, since (SA^) = 0, the minimal noise sum 
is uniform on the whole sphere 

^rnin ^-*max {^-^u) min 1 

The entanglement is maximal and constant in the (S2, S3) plane 
ls 2 ,s 3 = minX ajfe = (8Xl) min + 1 

a, b 

In this plane all modes have identical noise 



(5Xl{0)) = -{l + (5Xl{e))) 



(23) 
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This result is easily understandable if compared to the transmission of a squeezed 
beam by a beamsplitter. It is well known that part of the incident beam squeezing is 
lost to the reflected beam and leaks into the environment. 

We measured the noise spectra for all basis. We observed squeezing in the a + compo- 
nent, the a_ component being at the shot noise level. We checked that the spectra of 
the x, iy, +45, i(— 45) modes were all identical and squeezed by half the amount of the 
a + component squeezing, consistently with 



Although the x, y modes are both squeezed for orthogonal quadratures, the 
difference with the previous case is that they are now correlated. We verified this by 
measuring the entanglement in each basis. The results, displayed in Fig. show that 
the circular components are indeed uncorrelated, while the x, y modes are entangled, as 
well as the ±45 modes. 






Tin (s) 



Figure 12. Entanglement in different polarization basis, in the circular polarization 
case: (a) X x , y (0), (b) l a+tCr -(6), (c) I +45 _ 45 (6»). 
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4. Polarization entanglement 

Up to now, we have determined the correlation properties of the beam and shown 
that non-separable states, namely the ±45 modes in the high frequency limit, were 
produced in our system. We show in this Section that the quadrature entanglement 
demonstrated previously can be mapped into a polarization basis, thus achieving 
polarization entanglement. 



4-1. Definition and scheme 




Figure 13. Polarization-entangled beams generation set-up. Using the quarter-wave 
plates or not allows for measuring either Ss 2 or Ss 3 ■ 

Let us consider two spatially separated polarization modes a and (3. To each one 
we associate a set of Stokes parameters. One can extend the inseparability criterion for 
two modes to a pair of these operators [SJ EH U\ ■ Polarization entanglement is achieved 
when 

= IMS? + S p 2 f) + (6(S? + ^) 2 >] < I {[5?, Sfl)\ + \([SS, Sg\)\ (24) 

Because of the cyclical commutation relations between the Stokes operators, the 
criterion now depends on the polarization state of the beams (here, (Sf) and <Sf». In 
our case, we use two quadrature entangled modes a and b, mix them on a polarizing 
beamsplitter with an intense coherent beam B, polarized at 45° [Fig. H3] . The resulting 
beams, a and /3, are composed with modes A a and B y , and and B x respectively. The 
Stokes parameters are 



S? = A\A a - B\B y Sf = B\B X - A\A h 
S« = A a B\ + A\B y Si = A\B X + A b Bl 
S% = i{A a B\ - A{B y ) S§ = i{A\B x - A b Bl) 
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Denoting by a B} a b and 6 B the field amplitudes and the B-field phase, and 
assuming that the B-field is much more intense than the A-field (a^ > a a ,a b ), the 
two beams are orthogonally polarized: (Sf) = —(Sf) = —a B . The polarization 
entanglement condition (|24|) then reads 



Z*, < 2a% (25) 

On the other hand, the Stokes parameters fluctuations are proportional to the a, b 
modes quadratures 



5S« = a B 5X a (6 B ), 5S% = a B 5X b (6 B ) 

5S«= -a B SY a (6 B ), = a B 5Y b (6 B ) 

The inseparability criterion is thus directly related to the entanglement between 
the a, b modes 

Tip = a 2 B T ajb (6 B ) (26) 

If one locks the phase 9 B in order to obtain T a ^{6 B ) = mm e T ab (9) = T ab < 2, 
the beams are then polarization entangled. The S2 and S3 Stokes parameters can be 
measured using the right combination of plates and beamsplitter [TO] . 

4-2. Experimental results 

We have shown that the maximally entangled modes were the ±45 modes at the output 
of the cavity We thus insert after the cavity a half-wave plate to set A a = A +i5 and 
A b = A_ 45 . We then lock the phase 6 B and measure 



S S2 = ((SX x (9 B )f) , S S3 = {(SY y (9 B )) 2 ) (27) 
The results are shown in Fig. El and we observe 



1 

1 ••• o/3x2\ 



S S2 = —(6(S« + S%) 2 }^0.9Q<1 



Ss 3 = CT (^3 Q + ^T)^0.96<l 



so that 



2a B 



X s 

— if- = Ss 2 + Ss 3 = Z+4 5) j(_45) — 1-92 < 2 



This value is consistent with the quadrature entanglement measurement fl20D . 
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1 2 3 4 5 

Time (s) 

Figure 14. Normalized noises of S§ + S% (a) and £3 + S3 (b), when the phase 63 is 
locked. 

5. Conclusion 

Using the non-linearity of cold atoms, we have generated a quadrature-entangled 
beam. The maximally entangled polarization modes have been found following a 
general method to characterize entanglement in a two mode system. We have 
stressed the equivalence between our scheme and other continuous variable entanglement 
experiments in which the correlations are created via the mixing of independent beams. 
Moreover, a simple interpretation of the quantum properties of such a system was given 
in the Poincar sphere. To demonstrate the entanglement we have performed a direct 
measurement of the inseparability criterion j2j using two homodyne detections. We have 
then achieved polarization entanglement by mixing our quadrature entangled beam with 
an intense coherent field. Experimental evidence of this entanglement was given by the 
direct measurement of the Stokes operators noises for each beam. 

Appendix A. Uncorrelated basis existence 

Starting from a "correlated" polarization basis (a, b) ((5A a 5A b ) ^ 0), we prove the 
existence of the "uncorrelated" basis (u,v), such that 



{SAJA V ) = 



(A.l) 
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A u = cos ®A a - sin §e iuJ A b 
A v = sin §A a + cos $e ilJ A b 

with 

cos2 $= l ^ )|2 - | ^ 2)|2 S in2$ = « 
\M\' M ' AT 

w/iere M = (5 A 2 a ) (5 A a 5 A b )* + (5 A 2 b )* (5 A a 5 A b ) 
and M = ^/A\M\ 2 + [ \(5Al)\ 2 -\(5A 2 b )\ 2 } 2 
it is straightforward to see that the u, v modes thus defined satisfy (jA.lj) . 

Appendix B. Correlations in the Poincare sphere 

We give a brief demonstration of the general properties about the Poincare sphere 
enunciated in Sec. 12.21 Given the Poincare sphere defined from the uncorrelated basis 
(u,v), we calculate the entanglement and squeezing in different basis. We recall that 
the entanglement X ab between modes a and b only depends on the correlation term 
| (5A a 5A b )\, whereas the noise sum S ajb depends on the quantity | (5 A 2 ) \ 2 + \ (5 A 2 ) \ 2 . For 
a given polarization basis one has 

\{5A 2 a )\ 2 = p 4 (6A 2 u ) 2 + a 4 (6A 2 v ) 2 + 2a 2 f3 2 (5A 2 u )(5A 2 v )cos2<j ) (B.l) 
\(SA 2 b )\ 2 = a\5A 2 u ) 2 + (3*(6A 2 V ) 2 + 2a 2 f3 2 (6A 2 u )(6A 2 v ) C os2<p (B.2) 

For a and (3 fixed, |(<L4^}| and are maximal pour <fi = [tt] 

\(SA 2 a )\ max = p 2 (6A 2 u ) + a 2 (5A 2 v ) and \{5A 2 b )\ max = a 2 {5A 2 u ) + (3 2 {5A 2 V ) 
=► + \(SA 2 b )\), =0 = max [\(5A 2 a )\ + |(^)|] = (5A 2 U ) + (5A 2 V ) 

a,b 

The noise sum is thus minimal for all the linearly polarized modes with respect to 
u, v, i.e. in the plane (S[, S' 2 ) 



^si,s' — min^a.fc — (5X 2 ) min + (5X 2 ) min = S m ,,, 

a,b 



If the ellipticity increases (0^0), the noise sum increases, as can be seen from 
()B.1IIB.2|) . In the meridional plane (S' 2 , S%) the a, b modes satisfy (a = (3 = 1/V2) 

(m 2 a )\ 2 ) s >,s' 3 = (\(SA 2 b )\ 2 ) s ,, s , = 1 -{(5A 2 u ) 2 + (6A 2 v ) 2 + 2(6A 2 u )(6A 2 v )cosm} 
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This quantity is minimal for the circularly polarized modes a*,b* (0 = 7r/2): 



\(5A 2 a .)\ + 1(^)1 = \(SAl) - (6A*)\ = min{|(^)| + \(8A 2 b )\} 



a.b 



The noise sum then equals its maximal value S max . Assuming (5A^) > (5 A 2 ,), this 
value reads 



^a* ,b* (^A a *) m j n -|- (^dXfj,) min {8X u ) m i n -|- {SX v ) max — ^max 

More generally, one has 

^max = min{ (SX^min + (8X%) max , (5X 2 ) max + (§Xy) m i n } 

Let us now consider the entanglement. Optimal entanglement is obtained by 
construction for modes a*, b* 



1a*,b* — mmX aj b — (SX u ) m i n + (SX v ) m i n — S m j, 

a,b 



The entanglement decreases with the ellipticity. For the linearly polarized modes 
= 0), Eq. @ yields 

\(5AJA b )\ = a(3\(5A 2 u ) - (5A 2 V )\ 
X aj b reaches its maximal value for the u, v modes 

1 u ,v = -{(fiXDmin + (5X 2 ) max + (5X 2 ) min + (5X 2 ) max } = maxl^j 

Z a,b 



max i 



However, X a ^ is not constant in the equatorial plane; it is minimal (I c ,d = S 
for the 45°-polarized modes (a = (3 = l/v^2), denoted by c, d. 

Last, we would like to point out that optimizing the squeezing sum of two modes 
is not equivalent to optimizing the squeezing for one mode only. Finding the maximally 
squeezed mode is not trivial, since the u, v modes are not a priori independent. No 
condition holds on (5A U 5A\). The noise of one quadrature of mode a is given by 

(5X 2 (6)) = (3 2 (5X 2 (6)) + a 2 (5X 2 (6 - 0)> - 2a/3(6X u (0)6X v (6 - 0)> 
The correlation term can be written as 



(5X U (6)5X V (6 - 0)) = (5AJAl)e-^ + (5 A\5 A v ) e* 4 * 
= 2cos{(j)c-<l>)\{8A u 8Al)\ 
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with (pc the phase of (d~A u 5Al). The complex general solution takes a simple form 
if we assume that {5A U 5A\) is a real positive number (<pc = 0). The optimal value is 
then reached for = and one finally gets 



It is therefore possible to obtain a better squeezing on one mode than that of u, v. 
Note that, if (5A U 5A\) = 0, the u, v modes are independent and, as in our experiments, 
the best squeezing is that of one of the u, v modes. 
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